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Abstract 

The most tight conformation of the trefoil knot found by the SONO 
algorithm is presented. Structure of the set of its self-contact points is 
analyzed. 



1 Introduction 

Finding the best way of packing a tube within a box seems to be rather a gar- 
dening than a scientific problem. However, the optimal single helix, discovered 
in a computer simulation study of this problem, [jlj and S, proves to be ubiqui- 
tous in many proteins as their a-helical parts. It seems, as suggested in ||, that 
also the closely packed double helix appearing in the process of twisting two 
ropes together [Q have been already discovered and applied by nature. Labora- 
tory experiments allow one to observe in the real time how the optimal helices 
are formed in various systems e.g. the bacterial flagellas || or phospholipid 
membranes ||. 

Both processes, of packing the ropes and twisting them together, occur si- 
multaneously when a knot tied on a rope becomes tightened. The problem of 
finding the most tight, least rope consuming conformations of knots was inde- 
pendently posed and indicated as essential by different authors; for references 
see [Q. Knots in such optimal, most tight conformations are often called ideal, 
a term proposed by Simon || , and introduced into the literature by Stasiak Q . 
Ideal conformations minimize the value of the size-invariant variable A = L/ D, 
where L and D are, respectively, the length and the diameter of the perfect rope 
(defined below) on which the knot is tied. The only knot whose ideal confor- 
mation is known at present is the trivial knot (unknot). See Fig.|l|. Its length 
in the ideal, circular conformation equals ttD, thus A = tt. Finding the ideal 
conformation of a nontrivial knot is a nontrivial task. Initiated a few years ago 
search for the ideal conformations of nontrivial knots continues. 



Figure 1: Ideal unknot. 



One of the algorithms used in the search is SONO (Shrink-On-No-Overlaps) [fiol . 
SONO simulates a process in which the rope, on which a knot is tied, slowly 
shrinks. The rope is allowed to shrink only when no overlaps of the rope with 
itself are detected within the knot. When such overlaps occur, SONO modifies 
the knot conformation to remove them. If this is no more possible, the pro- 
cess ends. Unfortunately, ending of the tightening process does not mean that 
the ideal conformation of a given knot was found. The tightening process could 
have stopped also because a local minimum of the thickness energy was entered. 
The possibility that there exists a different, less rope consuming conformation, 
cannot be excluded. 

SONO has been used in the search of ideal conformations of both prime and 
composite knots. Parameters of the least rope consuming conformations found 
by the algorithm were listed in (ll[] and Jl^ . In a few cases, SONO managed to 
find better conformations than the simulated annealing procedure Q . However, 
for the most simple knots, in particular, the trefoil knot, the simulated annealing 
and SONO provided identical results; the A values are identical within experi- 
mental errors. It seems obvious, that no better conformations of the knot exist. 
We feel obliged to emphasize, however, that it is only an intuitively obvious 
conclusion - no formal proofs have been provided so far. As indicated in we 
are in a situation similar to that, which lasted in the problem of the best pack- 
ing of spheres for 400 years. That the face centered cubic and hexagonal close 
packed lattices were among the structures which minimize the volume occupied 
by closely packed hard spheres seemed to be obvious since the times of Kepler, 
however the formal proof of the conjecture was provided but a few years ago fu^f . 
Waiting for the formal proofs that what we have observed in the knot tighten- 
ing numerical experiments is the ideal conformation of the trefoil, seems to be 
a too cautious attitude. Thus, after a few years of experimenting, we decided 
to present the best, least rope consuming conformation of the trefoil knot we 
managed to find. We compare it with the most tight conformation of the knot 
which can be found within the analytically defined family of torus knots. In 
particular, we describe the qualitative change in the set of self-contacts which 
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Figure 2: The perfect rope. Perpendicular sections of the rope are of the disk 
shape. None of the disks are allowed to overlap. This puts a limit not only on 
the spacial distance of different fragments of the curve into which the rope is 
shaped, but also on its local curvature. 

takes place within the trefoil knot during the tightening process. We believe 
that some of the features of the self-contact set we have found may be present 
also in ideal conformations of other knot types. 

An alternative method of searching for the most tight conformations of knots 
consists in inflating the rope on which the knot has been tied In such a process 
the length of the rope is kept fixed. The maximum radius to which the rope 
in a given conformation of a knot can be inflated is closely related with the 
injectivity radius considered in detail by Rawdon fl4|| . 

2 The perfect rope 

It is the aim of the computer simulations we perform to simulate the tightening 
process of knots tied on the perfect rope: perfectly flexible, but at the same time 
perfectly hard in its circular cross-section. The surface of the perfect rope can 
be seen as the union of all circles centered on and perpendicular to the knot 
axis C. See Fig.|[ 

We assume that C is smooth and simple, i.e. self-avoiding, what guaranties 
that at each of its points r the tangent vectors r(r), and thus the circular 
cross-section, are well defined. The surface remains smooth as long as: 

A. the local curvature radius r K of the knot axis is nowhere smaller than 
D/2, 

B. the minimum distance of closest approach d* is nowhere smaller then D/2 

The minimum distance of closest approach d*, known also as the doubly 
critical self-distance, see is defined in |Tq| , as the smallest distance between 
all pairs of points (ri, Y2) on the knot axis, having the property, that the vector 
(v2 — i"i) joining them is orthogonal to the tangent vectors r(ri),r(r2) located 



3 



Figure 3: The trefoil knot is a torus knot - it can be tied on the surface of a 
torus. 

at the points: 

d*(C)= min {|r 2 -n| : r(n) 1_ (r 2 - n), r(r 2 ) i_ (r 2 - n)} (1) 

ri,r 2 SC 

As shown by Gonzalez and Maddocks |l6| , the two conditions can be gath- 
ered into a single one providing that the notion of the global curvature radius 
p G is introduced: 

PgOi) = min p(ri, r 2 , r 3 ) (2) 

r 2 ,r 3 6C 

where, p(ri,r 2 ,r3) is the radius of the unique circle (the circumcircle) which 
passes through all of the three points: ri,r 2 and r3. Using the notion of the 
global curvature, the condition which guaranties smoothness of the knot surface 
can be reformulated as follows: 

C. the global curvature radius p G of the knot axis is nowhere smaller than 
D/2. 

Analysis of the conformations produced by the SONO algorithm proves that 
conditions A and B, (and C) are fulfilled. 

3 Parametrically tied trefoil knot 

The trefoil knot can be tied on the surface of a torus. See Fig.||Consider the set 
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of 3 periodic functions: 



x = [-R + rcos(2^ 1 7ri)] sin{2v 2 irt) (3) 
y = [R + rcos(2v 1 irt)]cos(2v 2 Tit) (4) 
■ sin(2 v\ 7r t) (5) 



r s 



The trajectory determined by equations ||, f| and || becomes closed as t spans 
a unit interval. For the sake of simplicity we shall consider the [0, 1) interval. 
For all relatively prime integer values of v\ , v 2 equations |§[ |] and || define self- 
avoiding closed curves located on the surface of a torus. R denotes here the 
radius of the circle determining the central axis of the torus while r denotes 
the radius of its circular cross-sections. For the trefoil knot, frequencies v±, v 2 
equal 2 and 3, respectively. In what follows we consider knots tied on a rope; 
trajectories defined by equations ||, || and || determine position of its axis. 

The (1/1,1/2) and the (y 2 , v\) torus knots are ambient isotopic, i.e. they can 
be transformed one into another without cutting the rope on which they are tied 
[17 . As shown previously, the (2, 3) version of the trefoil is less rope consuming 
. Thus, the (3, 2) version will not be discussed below. 

Assume that the trefoil knot whose axis is defined by equations |[ || and || 
is tied on a rope of diameter D = 1. In what follows we shall refer to it as the 
parametrically tied trefoil (PTT) knot. In such a case, radius r of the torus on 
which the axis of knot is located, cannot be smaller than 1 /2 ; below this value 
overlaps of the rope with itself will certainly appear; at r — 1/2 the rope remains 
in a continuous self-contact along the torus axis. To keep the self-contacts we 
assume in what follows that r = 1/2. To check, if the knot is free of overlaps in 
other regions, one can analyze the map of its internal distances. Let t\ and t 2 be 
two values of the parameter t, both located in the [0, 1) interval. Let (x±,y±, z\) 
and (x 2 ,y 2 ,z 2 ) be the coordinates of two points indicated within the knot axis 
by t\ and t 2 , respectively. Let d(ti,t 2 ) be the Euclidean distance between the 
points: 



d(ti ,t 2 ) = VO2 - zi) 2 + (y 2 - yi) 2 + {z 2 - z x f (6) 

The map of the function, see Fig.|| displays a mirror symmetry induced by 
the equality d{t\, t 2 ) — d(t 2 , t{). 

Looking for possible overlaps within the knot one looks for regions within the 
internal distances landscape, where d{t\,t 2 ) < 1. The most visible depression 
within the landscape of the interknot distances is located around the diagonal 
where t\ — t 2 . As easy to see, d(t\,t 2 ) — along the line, but for obvious 
reasons this does not implies any overlaps within the knot. 

Another valley within which d(ti,t 2 ) may go down to the critical 1 value is 
localized in the vicinity of lines defined by equality \t 2 —ti\ = 1/2. To see, if in 
the vicinity of the lines the height really drops to or even below 1, we plotted 
the map of the d(ti, t 2 ) function in such a manner, that regions lying below the 
arbitrarily chosen 1.005 level were cut off. 
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Figure 4: The map of the intraknot distances of the most tight PPT knot. 



As seen in Fig|5j there are four such regions within the PTT knot: one in the 
shape of a sinusoidal band and three in shapes of almost circular patches. The 
band contains in its middle the mentioned above continuous line of self-contacts 
points; it is the axis of the torus on which the knot is tied. The circular patches 
contain 3 additional contact points; when R becomes too small, overlaps appear 
around the points. Numerical analysis we performed reveals that (with the 5 
decimal digits accuracy we applied) the overlaps occurring within these regions 
vanish above R = 1.1158. For R = 1.1159 the distance between the closest 
points located within these regions of the knot equals 0.9999. For R = 1.1158 
the distance is equal 1.0000. Where, within the PTT knot the self-contact points 
are located is shown in Fig.^ 

4 SONO tied trefoil knot 

Considerations presented above indicated the value of R, at which the PTT 
knot reaches its most tight conformation. The length L t of the rope engaged 
in this conformation of the trefoil knot equals 17.0883. Can one tie the trefoil 
knot using a shorter piece of the rope? Theoretical considerations indicate 
that this possibility cannot be excluded. As proven in the piece of rope 
used to tie the trefoil knot cannot be shorter than Lm = (2 + Vfyn ~ 10.72. 
Such a location of this lower limit leaves a lot of place for a possible further 
tightening of the knot. Application of SONO reveals that the tightening is 
possible providing the conformation of the knot is allowed to leave the subspace 
of the parametrically tied torus conformations. This happens spontaneously in 
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Figure 5: The map of the intraknot distances. Left - the most tight PTT knot. 
Right - the most tight STT knot. The map was cut from below at the height 
10.005. 



numerical simulations in which the most tight PTT knot is supplied to SONO as 
the initial conformation. SONO algorithm manages to make it shorter. In the 
simulations we performed, SONO reduced the length of the knot by about 4% 
to L exp — 16.38. The discrete representation of the knot used in the simulations 
contained N — 327 nodes. Below we describe the final conformation. For the 
sake of simplicity we shall refer to trefoil knots processed by the SONO algorithm 
as the SONO tied trefoil (STT) knots. 

The differences in the conformation of the most tight conformations of the 
PTT and STT knots is a subtle one. The essential difference lies in the structure 
of the sets of their self-contact points. As mentioned above, the circular line of 
self-contact points present in the family of the PTT knots stays intact as R is 
changed within the family. Tightening of a PTT knot achieved by decreasing the 
radius R of the torus stops when additional discrete points of contacts appear 
at three locations within the knot. This happens as R becomes equal 1.1158. 
Further tightening of the knot within the family of PTT knots is not possible, 
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Figure 6: Localization of the set of self-contact points within the most tight 
PPT knot. 



it becomes possible within the family of the STT knots. 

During the tightening process carried out by SONO, the set of the self- 
contact points undergoes both qualitative and quantitative changes. First of 
all, the line of contacts present in the PTT knot changes its shape becoming 
distinctly non-circular. Secondly, the three contact points give birth to pieces 
of new line of self-contacts. Unexpectedly, the new pieces do not connect into a 
new line, wiggling around and crossing the old line, but they are mounted into 
the old line in such a manner, that a single, self-avoiding and knotted line of 
self-contacts is created. That this is the case was revealed by a precise analysis 
of the intcrknot distances function. A map covering the interknot distances only 
within the very thin [1.00000, 1.00002] interval shows two separated lines, see 
Fig.^, corresponding to a single, self-avoiding and knotted line of contact. 

In addition to the line, a set of three points of self-contacts is formed. The 
points are located at places where the line of self-contacts becomes almost tan- 
gent to itself. The self-contact line runs twice around the knot. As a result, 
each of the circular cross-sections of the rope stays here in touch with another 
two such sections. The close packed structure formed in such a manner is much 
more stable than the structure of the most tight PTT knot, where single con- 
tacts were predominant. Let us note, that figure le presented in ref. M] a 
similar self-contact line structure can be seen. Unfortunately, inspecting the 
figure one cannot see, if the "self-contact spikes" shown there form a single, 
self-avoiding, knotted or a double, crossing itself line. The problem was not 
discussed in the text. Let us emphasize, however, that the difference between 
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Figure 7: The set of the self-contact points in the most tight STT knot as seen 
within the map of the intraknot distances. 

the two possibilities is confined to a zero-measure set. 

5 Discussion 

Ideal knots are objects of which very little is known still. The only knot whose 
ideal conformation is known rigorously is the unknot. Its ideal conformation, a 
circle of a radius identical with the radius of the rope on which it is tied, can 
be conveniently described parametrically. The set of the self-contact points is 
here limited to a single point: the center of the circle. All circular sections of 
the rope meet at this point. The maximum local curvature and the minimum 
double critical self-distance limiting conditions arc simultaneously met. 

The situation in the case of the trefoil knot, the simplest non-trivial prime 
knot, is radically different. Here the most tight parametrically defined con- 
formation proves to be not ideal. As demonstrated by the present authors, it 
can be tightened more with the use of the SONO algorithm. The set of the 
self-contact points becomes rebuilt during the tightening process. Its topology 
becomes different. In the case the PPT knot the set of the self-contact points 
consists of acircle and 3 separated points. As the numerical experiments we per- 
formed suggest, in the case of the STT knot, the set of the self-contact points 
turns unexpectedly into a single line. Which the structure of the set of self 



9 



Figure 8: Position of the line of the self-contact points within the ideal trefoil 
knot. To make the line more visible, a part of the knot was cut out. 



contact points in other prime knots is, remains an open question. 
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